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We have studied the pseudogap formation in the single-particle spectra of the half-filling two- 
dimensional Hubbard model. Using a Green's function with the one-loop self-energy correction of the 
spin and charge fluctuations, we have numerically calculated the self-energy, the spectral function, 
and the density of states in the weak-coupling regime at finite temperature. Pseudogap formations 
have been observed in both the density of states and the spectral function at the Fermi level. 
The pseudogap in the spectral function is explained by the non-Fermi-liquid-like nature of the self- 
energy. The anomalous behavior in the self-energy is caused by both the strong antiferromagnetic 
spin fluctuation and the nesting condition on the non-interacting Fermi surface. In the present 
approximation, we find a logarithmic singularity in the integrand of the imaginary part of the self- 
energy. Concerning the energy dependence of the spectral function and the self-energy, two theorems 
are proved. They give a necessary condition in the self-energy to produce the pseudogap at the Fermi 
level. The pseudogap in the spectral function is highly momentum dependent on the Fermi surface. 
It opens initially in the (±7r, 0), (0, ±7r) regions as the normal state pseudogap observed in the high- 
T c superconductors and if the interaction is increased, it spreads to other Fermi surface sectors. 
The anisotropy of the pseudogap is produced by the low-energy enhancement of the spin excitation 
around Q = (7T,vr) and the flatness of the band dispersion around the saddle point. 

PACS numbers: 71.10.Fd,71.20.-b 



I. INTRODUCTION 



The recent series of angle-resolved-photoemission experiments in high-T c superconductors have brought a large 
number of new insights concerning the low energy single-particle states of these materials. The experimental ob- 
servations include a rather wide flat-band region around (tt, 0) in the effective band dispersion and also the 
anisotropy of the normal state pseudogap which is consistent with the d-wave symmetry superconducting gap. 
Many theoretical approaches and computer simulations have been performed to investigate the pseudogap. Naturally, 
our ultimate objective is to explain the existing experimental data. However, the main purpose of this work is not 
to give a direct explanation of the experiments, but rather, to study the pseudogap formation in the Hubbard model 
which is often used to describe the high-T c materials. 

To avoid confusion, we should make clear what kinds of pseudogap we are dealing with in this paper. In reality, in 
the literature this term has been used to represent several different features. What it is common in all cases is the 
fact that the pseudogap indicates a suppression or disappearance of the spectral intensity at the Fermi level. We can 
list the following three different cases. The first one is a disappearance of the spectral intensity due to the downward 
shift of the band dispersion around (tt, 0) region Jl(J. For example, some of the quantum Monte Carlo (QMC) 
simulations show such behavior in the strong coupling regime and at very low dopings |^d|-^3|. This phenomenon 
is quite interesting since it looks as if it violates the Luttinger theorem on the Fermi surface. The second case is 
an effective pseudogap in the spectral function between a strong quasiparticle peak and a weak satellite peak. It is 
possible to find this situation even in a normal Fermi liquid. The third one is a suppression of the single-particle 
peak at the Fermi level with a two peak structure in its place. The pseudogap of our interest is this latter case. It 
should be clearly distinguishable from the second case. One physical origin of this kind of pseudogap is of course the 
fact that it is a precursor of the spin density wave (SDW) transition which takes place in the 2-d Hubbard model for 
zero doping at T = 0. We are aware of some models which take into account precursor effects of the superconducting 
transition to explain the pseudogap behavior observed in the normal phase of high-T c materials |p^-|l7t . We don't 
consider the superconducting transition effects in this paper. 

In our work, we provide an explanation of the origin of the pseudogap formation in the 2-d Hubbard model, at 
half-filling, in the weak coupling regime. To calculate the spectrum we use a single-particle Green's function which 
has a paramagnon-like one-loop self-energy correction for both charge and spin channels. At first we show the results 
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of our numerical calculations of the density of states (DOS), the spectral function, and the self-energy. The analysis 
of the real and imaginary parts of the self-energy plays a crucial role in the understanding of the microscopic origin 
of the possible structures of the spectra. Following this, we discuss the detailed origin of the pseudogap formation 
within our formulation. 

We found pseudogap formations at the Fermi level in both the spectral function and the DOS. The pseudogap in the 
spectral function is explained by the non- Fermi- liquid behavior of the self-energy, i.e., its imaginary part has a negative 
peak and its real part has a positive slope. We show that both the strong antiferromagnetic spin fluctuation and the 
nesting condition on the non-interacting Fermi surface are the origins of the anomalous behaviors in the self-energy. 
We also generally argue the relation between the spectral function and the self-energy at the Fermi level. We show, 
for example, that if the spectral function has a pseudogap, the real part of the self-energy has a positive slope at the 
Fermi level. An auxiliary relation is used to present the argument. This relation derived from the Kramers-Kronig 
transformation is proved in the appendix. 

The numerically obtained pseudogap shows a strong momentum dependence. For a certain choice of the parameters, 
the pseudogap appears around (7r, 0) but not at (tt/2,it/2). We show that the anisotropic pseudogap comes from the 
low-energy enhancement of the spin excitation around Q = (tt, 7r) and the flatness of the band dispersion. 



II. MODEL 



In the momentum representation the Hubbard model, with a nearest-neighbor hopping, on a two-dimensional square 
lattice can be written as 



U 



p 2 -k-0' a P2-0-apia' 



(2.1) 



k,pi,p 2 ,(T 



with £k = — 2t(cosk x + cosfcj,), a representing the spin index and with U being the coupling constant of the on-site 
Coulomb repulsion between two electrons with opposite spins. We redefine the chemical potential fi by adding the 
Hartree-Fock term and this is chosen as the origin of the energy for the half-filling case. 

We next introduce the single-particle Green's function at finite temperature in the paramagnetic regime. In general, 
we can write the full single-particle Green's function as 



G(k,^„)- 1 = Go(k, w „)~ 1 - E(k,w n ) 



(2.2) 



where Go(k, u n )~ 



£k and E(k, u> n ) is the proper self-energy. Here, in the weak-coupling regime, we con- 



sider the one-loop self-energy correction to the Green's function. The formulation we use is based on the standard 
perturbation expansion over U/t. These one-loop corrections are constructed with the free (mean-field) fermion line 
and the corresponding charge or spin fluctuation lines. The charge and spin fluctuations are described in terms of 
the susceptibilities calculated within the random phase approximation (RPA). We neglect higher-order contributions 
which appear in the diagrammatic series expansion. This kind of approximation has already been used in paramagnon 
theory pl|. Thus, the one- loop self-energy is simply 



£(k,c n ) = [/ 2 T]T£ 



:Xc(q, v m ) 



+ 2*s(q, fm) -xo(q,^m) 



G (k + q,o; n +v m ) 



(2.3) 



where we have defined the charge susceptibility as X c (q, v m ) = Xo(q> v m)/\\ + ^Xo(q> v m)] ancl the spin susceptibility 
as Xs(q, v m) = Xo(q> v m )j\Y — ?7xo(q, v m )). To avoid the double counting of the second-order self-energy diagram, we 
have subtracted xo(q> v m) in the square brackets in E(k, u> n ). Here xo is given by 



v(n „ ^ _ y /(£k+q) - /(£k) 

A -fiv m - (e k +q - £k) 



(2.4) 



with v m — 2mrntT and f(x) is the Fermi distribution function defined by f(x) = l/[exp(x/T) + 1]. 
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III. NUMERICAL RESULTS 



To calculate the spectra obtained from the above Green's function, we made the analytic continuation of the energy 
variables of the self-energy and the Green's function. In this way we obtained the retarded self-energy and the retarded 
Green's function. As a result, the imaginary part of the self-energy becomes, 



ImS fl (k,cj) = -U 2 J2 



-Imxf (q, £ k+q 
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-Imxc (q, £k+ q -u})- Imxo (q, e k+q - lo) 



X [/fck+q) + «(£k+q - u)] 



(3.1) 



where n(x) is the Bose distribution function defined as n(x) — l/[exp(x/T) — 1]. 'R' denotes the retarded function. 
The corresponding real part can be obtained from the imaginary part by means of the Kramers-Kronig transformation. 
Using those schemes, we have numerically calculated the susceptibilities, the self-energy, and the Green's function. 
The momentum summation in Imxo(qi v ) wa s reduced to a contour integral in the Brillouin zone and roughly 4000 
line-elements have been numerically summed up along the contour. The momentum summation in the imaginary part 
of the self-energy and the DOS have been done on a 120 x 120 mesh in the Brillouin zone. All parameters used in the 
calculation are within the Stoner instability condition in RPA. 



A. Density of States 

In Fig. |l| we show the U dependence of the DOS N(co) — ^2 k A(\t,u]) for T/t — 0.5 where the spectral function 
A(k, lo) is defined by 

A(k,w) = --lmG R Ck,u) 

1 ImE' R (k,w) 
~ ~n [lo - e k - ReS fl (k, lo)} 2 + [ImS«(k, lu)} 2 ' ^' 2 > 

For U/t — 1.5, the overall structure still keeps the shape similar to the non-interacting case with a clear single peak 
on the Fermi level. For U/t = 2.5, a weak pseudogap can be seen, and finally, for U/t = 3.0, a rather wide pseudogap 
can be observed. The high-energy band tails are developed as U/t becomes large. 

Figure |2| shows the temperature dependence of the pseudogap formation in the density of states for U/t — 2.0. The 
pseudogap develops as we lower the temperature. Using the Stoner condition, we numerically found that T/t = 0.21 
is the SDW transition temperature for U/t = 2.0. The temperature effect appears only around the Fermi level and 
the overall high-energy structure is left almost unchanged. 

A similar calculation of the DOS was done in the work by Kampf and Schrieffer jl9j . They have also obtained the 
pseudogap formation near the Fermi level lo = for a weakly doped Hubbard model at zero temperature. In their 
calculation they took into account only the spin fluctuation effects in an one-loop self-energy similar to ours. Since the 
spin fluctuation effect becomes much stronger than the charge contribution near the SDW transition, the similarity of 
our calculation and theirs is natural. Pseudogap openings in the DOS in the weak-coupling regime has been observed 
in a QMC simulation [2(J. They obtained pseudogaps for U/t = 4.0 and T/t= 0.2 for several lattice sizes from 8x8 
to 12 x 12. A similar behavior in the DOS obtained by the fluctuation exchange (FLEX) approximation They 
obtained a very narrow and deep pseudogap in the DOS for U/t = 1.57 and T/t = 0.05 at half-filling. All of those 
results including ours are at least qualitatively in good agreement with each other. The DOS for U/t = 3.0 in Fig. |l| 
has a rather large pseudogap. This may be due however to our simple approximation. 



B. Spectral function and the self-energy 

The existence of the pseudogap in the DOS enables us to guess the decrease of the intensity of the spectral function 
at lo = 0. In Fig. ^j, we plot the temperature dependence of the spectral function A(kp,Lo) at hp — (n,0) together 
with the real and imaginary parts of the self-energy. We have also calculated A(k, lo) for other momenta and our 
results are a strong indication of the anisotropy of the pseudogap. We will discuss it in more detail in a later section. 

As we have expected, Fig. || (a) shows the weight decrease in the spectral function as the pseudogap develops in the 
DOS shown in Fig. ||. The self-energy shows a specific behavior. The real part (Fig. || (b)) has a positive slope and the 
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imaginary part (Fig. [| (c)) shows a negative peak. This tendency becomes stronger as we decrease the temperature. 
Moreover, the Green's function has three poles which correspond to the solutions of u> — £k F = KeY, R (kp , u>) for 
T/t = 0.22 as shown in Fig. || (b). Two of the poles are associated with two peaks of A(kp,u>) in Fig. || (a) and 
the other pole is linked to the pseudogap. Those features in the spectral function and the self-energy clearly indicate 
the destruction of the Fermi-liquid quasiparticle states. We will argue in more detail on the relation between the 
pseudogap and the anomalous behavior in the self-energy in the next section. 

This same trend in the spectral function has already been shown in the so-called two-particle self-consistent (TPSC) 
approach by Vilk and Tremblay p2|p3| . In the FLEX approach |2l), however, there is no pseudogap in the spectral 
function although their self-energy shows a similar non-Fermi-liquid behavior that both our result and TPSC approach 
have shown. The detailed analysis and comparisons between those approaches can be seen in Ref. ]23| . 

For the weak coupling regime at half-filling, the finite temperature pseudogap formation in A(k, u) in a QMC 
simulation pij| was attributed to a finite size effect. Making a different analysis of the QMC simulation data, the 
pseudogap formation was confirmed by Creffield et al pp|] . The latter authors have used a singular- value-decomposition 
method instead of the maximum-entropy method to obtain the spectral function from the calculated finite temperature 
Green's function. They obtained a clear pseudogap opening even with a 12 x 12 system size. 



IV. PSEUDOGAP AND THE NON-FERMI-LIQUID BEHAVIOR IN THE SELF-ENERGY 

In this section we discuss the origin of the non-Fermi-liquid behavior in the imaginary part of the self-energy. This 
section is organized in three subsections. Firstly we analyze the imaginary part of the self-energy within the present 
paramagnon theory. A log singularity is found in the integrand of the imaginary part of the self-energy. Next, we 
examine a possibility of having the same tendency of the imaginary part of the self-energy using a model susceptibility 
which does not produce the log singularity. Finally, we discuss a general relation between the self-energy and the 
spectral function around the Fermi level. The argument in the final subsection does not depend on the models used 
and on the origin of the pseudogap. 



A. Log singularity in the integrand of the imaginary part of the self-energy 

We can now explain the origin of the non-Fermi-liquid negative peak behavior in Im£(k, oj) in the present formu- 
lation. The point is how to take into account the strong enhancement of the spin fluctuation around q = (n, tt) = Q. 
Let's consider the integrand of the spin component of the self-energy at k = kp. The integrand is 

3 

4 F (q,w) = --U 2 lm X R (q,£ kF +q - [/(skj,+q) + n(£ kj ,+ q - w)] . (4.1) 

Considering the strong low-energy enhancement of Imx s (Q, v), we suppose that the major contribution to the 
imaginary part of the self-energy comes from the integrand at q = Q. By setting q = Q, the spin fluctuation 
spectrum can be written as 

Imxf (Q, £kF+Q - «) = ImX " fl( 2 Q '" ) j (4.2) 

[l-t/Rcxo fi (Q,--)] 2 +[t/Im Xo «(Q,^)] 2 

with 

ImxflQ, v) = |po (|) tank (£) (4-3) 

R erf<Q,„> = ^Il4<2^. ,4.4) 

Here po{x) — (l/(27r 2 t))K (y/ 1 — \x] (4t)] 2 ^ is the density of states for the non-interacting band £k where K(fc) is the 
complete elliptic integral of the first kind. In the same way, we also have 

/(£k F +q) + ™(£k F +q - w) = -- coth (|^) . (4.5) 



and 
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Note that we have used the nesting condition £k F +Q = — £k F = on the Fermi surface. Then, we obtain 

I (Q uj) = 3 L ,2 -W4)poW2)tanh(^/(4T))coth^/(2r)) 
" 2 [l-C/Rcx fl (Q,-^] 2 +[t/Im X «(Q,- W )] 2 ■ 

For uj ~ 0, by the cancellation of the cj/T dependence in tanh and coth, we see that the numerator always has a 
log singularity since 

_ l po (0L) tanh ( Z0L) coth ( JL) „ JL In (™ V (4.7) 
4 P0 \2/ V 4 ^; V2T/ 167rf V M / 

By approximating the denominator with its value at w = 0, the behavior of the integrand for q = Q is dominated by 
the log singularity, 

^(Q,*)—|£ ln(32VN) a - (4-8) 

In lowering the temperature or by increasing U, 1 — t/Rcx^(Q,0) becomes smaller and the contribution from the 
region around Q = (tt,tt) dominates the other region's contribution and it causes the larger peak of ImE(ki?,cj) at 
uj = 0. 

In Fig. H we plot the temperature dependence of I^ f (Q,uj) for U/t — 2.0 using Eq. ( f4.6| ). Fig. || (a) shows the 
low-energy enhancement of Imx s (Q,i / ). The log divergence at uj — exists for any finite values of T. /k F (Q,w) 
shows a remarkable enhancement at the same time as Imx s (Q,^) has a sharp increase. Clearly, one can see that 
this enhancement of ik F (Q,w) around uj = causes the non- Fermi-liquid, i.e., the negative peak structure in the 
imaginary part of the self-energy. 



B. No log-singularity case 

In the last subsection we have found a logarithmic behavior in the integrand of the imaginary part of the self-energy. 
This has a strong influence in producing the non-Fermi-liquid behavior in the self-energy. In general, however, the 
logarithmic van Hove singularity in the DOS will be smeared to some extent by any finite interaction between 
electrons. One might suppose that the non-Fermi-liquid-like behavior we have seen in the last section is an artifact of 
our approximation. In this subsection, we consider this question by introducing a model susceptibility which does not 
possess the log-singularity behavior. Again we take into account only the spin components and neglect other channels 
in the self-energy. 

From Imx^(Q, v) given by Eq. ( |4.3| ), we see that the logarithmic behavior comes from the non-interacting density 
of states po(ui). To consider the case that Po(uj) does not have any divergence at uj = 0, we simply assume that po(ui) 
is finite at uj = 0. Then, from Eq. ( |4.3| ), we can approximate the imaginary part of Xo(Qj^) as Imxo^Q,^) = cv 
around v = where c is a positive constant. Since we would like to sec the low-energy behavior of the integrand of 
ImE, it is sufficient to approximate the spin susceptibility around v = 0. Furthermore, approximating Rexo(Q: v) by 
its value at v = 0, we obtain the imaginary part of the model susceptibility written as 

Im X ; fl (Q, v) = 5 . (4.9) 

[l-E/Re X £(Q,0)] 2 + M2 

One can see that our model susceptibility is formally equivalent to the phenomenological model susceptibility in 
the nearly-antiferromagnetic- Fermi-liquid (NAFL) theory |25|-|28|]. This model was introduced by Millis, Monien, and 
Pines to explain the NMR experiments in the high-T c materials [^9| . By applying the NAFL theory to the effective 
models including the next-nearest-neighbor hopping, anisotropic pseudogap formation p7[ and band dispersion |2§| ] 
which explains the experimental data in the high-T c materials have been obtained. 

Using our model susceptibility, the integrand of the imaginary part of the self-energy becomes 

OQ,") = -fo ( C /2Wh(,/( 2 T)) 

2 [l-t/Rex£(Q,0)] 2 + M 2 

In Fig. H (a), we show the U dependence of ImXs R (Q, v) for T/t — 0.25. The necessary parameters in the model 
susceptibility are chosen as c/t — 0.5 and Rexo(Q,0) = 0.4. We have roughly evaluated those parameters from the 
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numerical data, but the accuracy is not important here. As U/t becomes large, Imx* ii (Q, v) shows a low-energy 
Stoner enhancement. In Fig. |5| (b) and (c), we show the U evolution of i£ (Q,a;). For small U/t, we see the u 
dependence similar to the Fermi-liquid-like lo 2 behavior around u> = although J£ (Q,w) is the integrand of the 
imaginary part of the self-energy. As U/t increases, this behavior disappears and finally a very sharp negative peak 
structure develops at to = 0. 

Thus, even if the density of states does not have a log singularity, the strong negative peak appears in the integrand 
of the imaginary part of the self-energy as the parameter set approaches the Stoner instability condition. 



C. Necessary condition for the pseudogap formation 

As our numerical results have shown, the non-Fermi-liquid behavior of the self-energy and the pseudogap are related 
to each other. In all cases, as we have seen in our results, a pseudogap in A(k.p,uj) accompanies the positive slope 
of the real part of the self-energy. In this subsection, we generally argue the relation between the self-energy and 
the spectral function at the Fermi level. We prove two theorems which hold between the self-energy and the spectral 
function. These theorems determine which conditions in the self-energy are necessary in order to have a pseudogap 
in the spectral function. Through the argument, we assume that: (i) the imaginary part of the (retarded) self-energy 
is always negative and finite, (ii) at u = 0, we have (£k F — m) + ReE^kp, 0) = 0. We leave the chemical potential [i 
to keep the generality. Those assumptions are physically reasonable. 

Next, we enunciate the two theorems, after that, we give the proofs of them. 

Theorem I — If the imaginary part of the self-energy has a maximum at the Fermi level, the spectral function has 
a maximum at the Fermi level. 

One easily sees that the conventional Fermi liquid satisfies the theorem I. 

Theorem II — If the spectral function has a minimum at the Fermi level, the real part of the self-energy has a 
positive slope at the Fermi level. 

Note that if the spectral function has a hollow at the Fermi level, we can say that the spectral function has a two 
peak structure around the hollow since the spectral sum over energy is finite. Hence, to discuss the existence of the 
pseudogap, it is sufficient to observe if the spectral function has a maximum or minimum at the Fermi level. 

Proof of theorem I — If ui is slightly different from 0, from the definition of the spectral function, we have 

A(k F ,u>) = - 3 . (4.11) 

TT |Im£«(k F ,w)| + [u - (e kF -/x) -Rel2 R (k F ,uj)f /\lm^ R (k F ,uj)\ 

If u) = Q, by the assumption (ii), we obtain 

A ^ = A^«^M (4 ' 12) 

From the condition of this theorem, Im£^(kp,w) has a maximum at the Fermi level, and we have |Im£ (kp, cj)| > 
llmE^kp, 0)|. Moreover, the second term of the denominator of Eq. (4.11) is a positive finite, then we find A(k F , 0) > 
A(k F ,oj). The theorem I has been proved. 

Before going to the proof of theorem II, we give the following lemma which is used in the proof. 

Lemma — If ImH R (k F , lo) has a maximum (minimum) at the Fermi level, then ReS^kp, ui) has a negative 
(positive) slope at the Fermi level. This lemma can be proved using a property of the Kramers-Kronig relation 
between lmI] R (k F , lo) and ReE fl (kp, oj). We give a proof of this property in the appendix, and the lemma follows 
immediately from it. 

Proof of theorem II — From Dyson's equation, we obtain that 

Re£*(k, ui) = w - (e k - fi) - ^G R (k,oj) _ ^ 



and 



[ReG R {k,oj)Y + [lmG R (k,oj)Y 



Im^(k,c) = Im f (k '" } - 2 . (4.14) 

[ReG R ik,uj)} 2 + [ImG R (k,Lj)f 



At k = kp, the latter equation gives 
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\lmE R (k F ,u;)\ = 5 . (4.15) 

\lmG R (k F ,oj)\ + [ReG R (k F ,uj)f /\lmG R (k F ,uj)\ 



At oj = 0, from the assumption (ii) and Eq. (4.13), we find that KeG R (k F ,0) = 0. Then, from Eq. (4.15), we have 



|I m ^(k r ,0)|= |ImGfl(k ^ o)| . (4.16) 

Since the spectral function has a minimum at oj = 0, we see \lmG R (k F , 0)| < llmG^ki?, cj)| for oj near zero. Thus, 
we have 

ImE^k^O) < lmS R (k F ,cj) < 0, (4.17) 

i.e., Im£, R (kp , uj) has a minimum at uj = 0. From the lemma, we see that Re'S R (k F ,uj) has a positive slope at oj = 0. 
The theorem II is finally proved. 

We easily see that the reverses of the two theorems do not always hold. For example, when ReY, R (k F ,uj) has 
a positive slope, A(k F ,uj) may have a single peak at the Fermi level. However, we can never have a simultaneous 
maximum in lmSi R (k F , uj) and a pseudogap in A(k F ,uj) at the Fermi level. The argument we gave here does hold 
very generally and it does not depend on the models and the origin of the pseudogap. 

The reverse of the lemma is not always true either. However, if we include an additional condition in the initial 
assumption such as the imaginary part of the self-energy has always a minimum or a maximum at the Fermi level, the 
revere of the lemma is satisfied. To prove this let's assume this new condition. From the lemma, ReY.(k F , uj) always 
has a (positive or negative) slope at the Fermi level. Thus, the maximum (minimum) of lmY. R (k F , oj) is automatically 
linked to the negative (positive) slope of ReY; R (k F , uj). The reverse of the lemma follows from this. 

By combining the reverse of the lemma and the theorem I, we immediately find that when ReYj R (k F ,uj) has a 
negative slope at uj = 0, the spectral function has a single peak at the Fermi level. In other words, as far as the 
Fermi-liquid-like negative slope remains in KeY] R (kp,uj) at oj = 0, the pseudogap is never produced in the spectral 
function at the Fermi level. 



V. BAND DISPERSION AND ANISOTROPIC PSEUDOGAP 

In Fig. U (a) we show the band dispersion extracted from the numerical data of the spectral function ^4(k, uj) for 
several momentum directions. The dot indicates the peak of the spectral function. The peaks roughly follow the 
original non-interacting band dispersion except around (ir, 0). Clear (pseudo)gap structure can be seen around (tt, 0). 
In contrast, around (7r/2, it/2), there is no such pseudogap structure. This exemplifies the anisotropy of the pseudogap 
on the Fermi surface. We show the structures of the spectral functions around (it, 0) in Fig. ^ (b). We chose k = (it, 0), 
(0.97r,0), and (0.857r,0) along the k x axis. By varying k, the symmetric two peak structure at k = (tt,0) gradually 
changes into a single peak structure with a satellite peak. The satellite peak corresponds to the shadow band of the 
SDW band dispersion. In our earlier work pOfl , another sign of similar shadow band formation around k = (tt/2, it/2) 
and (tt, 0) has been observed in the dispersion spectrum for a weakly doped Hubbard model. 

In Fig. we compare the spectral functions and the self-energies for kp — (tt,0) and (tt/2, tt/2). ReT, R (k F ,uj) at 
kp = (tt, 0) in Fig. (b) has a large positive slope around oj = 0. ReS fl (ki?, oj) at kp — (tt/2, tt/2) also has a positive 
slope, but it is still weak. lmE R (kp, uj) at kp — (tt, 0) in Fig. [j] (c) shows a sharp negative peak at oj = 0. In contrast, 
lmS R (kp, uj) at kp = (tt/2, tt/2) has a weak minimum at uj = 0. These remarkable momentum dependences of the 
real and imaginary parts around oj = determines the opening of the pseudogap. 

For larger values of U/t, we have obtained a pseudogap in A(k,uj) also at kp = (tt/2, tt/2). In Fig. |[ we show the 
spectral function for U/t = 3.0 and at T/t = 0.5. We chose the higher temperature to avoid the Stoner instability. 
The pseudogap at (tt/2, tt/2) is weaker than the one at (tt, 0). Comparing with the spectral function for U/t = 2.0 
shown in Fig. ^ (a), the intensity of A(kp,uj) in Fig. ^ is weak and its width are rather wide. This indicates band 
broadening effects due to the larger value of U/t. 

To understand the reason for the anisotropy of the pseudogap, it is sufficient to analyze the imaginary parts of the 
self-energy at oj — for kp = (tt, 0) and (tt/2, tt/2). There are two factors: the difference of the effective region of the 
q summation and the different behaviors of the bosonic excitation energy in the imaginary part of the self-energy. The 
first one is a kind of the selection rule of the momenta of the fermion-boson excitation and the latter is the restriction 
of energy on its excitation. We explain them as follows. 
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The distribution function term /(£k F + q ) + n(£k F +q — produces an important restriction to the q summation in 
Im£(kF, 0). In Fig. |^ the approximated summation area in the (q Xl gj,)-plane is shown by the shaded region in which 
q satisfies the condition |/(£k F + q ) + «(£k F + q )| > 0.5 for T/t = 0.22. We see that for both kp — (tt,0) and (tt/2, tt/2) 
the region around q — Q contributes to the q summation. However, the shaded area around Q for hp — (7r,0) is 
greater than that for (7r/2,7r/2). Thus, ImE(k^ ,u) for hp — (tt, 0) has a larger contribution from the spin excitation 
around Q = (tt, tt). 

The bosonic energy £ b (q) = £k F + q - w in the expression of lmS(kp, ui) for hp = (tt,0) behaves in a totally 
different way from that for hp = (tt/2, tt/2). To see the behavior for co ~ and around q ~ Q, we can write 
£ b (Q + Sq) ~ — £k F +5q, where <5q is a small vector to represents deviation from Q. Thus, we see that the behavior of 
£ h (Q + <5q) is determined by the band dispersion £k around k = kp. Hence, the bosonic energy contribution to ImE for 
hp = (7r,0) is quite concentrated within a small region around energy near the saddle point (flat-band) dispersion 
near (tt, 0). In contrast, the bosonic energy for hp = (ir/2, tt/2) can easily deviate from as q moves in the summation 
region near Q since £k F has a linear slope around hp = (7r/2,7r/2). Thus, the low energy enhancement in the spin 
fluctuation around Q makes a stronger contribution to lmE,(hp,uj) at hp = (ir,0) than that at hp = (tt/2, tt/2). As 
a result, the self-energy has a remarkable momentum dependence, and the anisotropy of the pseudogap is explained. 

There are some QMC simulations of the spectral functions in the weak-coupling regime at half-filling |2(||24|,|3l| . 
With the exception of the results from Creffield et al. (Ref. J2(|), a clear evidence of the anisotropic pseudogap in the 
spectral function is lacking in the QMC simulations. We think that this is also due to the cluster size effect and the 
difference in method in extracting the information of the spectral function from the Green's function as emphasized 
in Ref. ^0|. We know of no other example of the anisotropic pseudogap in the t-U Hubbard model. However, in 
FLEX calculations, for the doped Hubbard model [ j32|]33[| , a clear momentum dependence in the imaginary part of the 
self-energy at to = near the Fermi momentum was obtained. They found that |ImE(ki?,0)| has a large value along 
the Fermi surface which becomes larger and larger as hp approaches the (tt, 0) region. This tendency is consistent 
with our results at half-filling shown in Fig. [?] (c). 

VI. CONCLUDING REMARKS 

In this work we have studied the pseudogap formation in the two-dimensional Hubbard model at half-filling. We 
obtain a pseudogap formation in the density of states and in the spectral function. The pseudogap of the spectral 
function is produced when two conditions are satisfied: (i) a strong anti-ferromagnetic spin fluctuation, (ii) a nesting 
condition on the Fermi surface. The pseudogap on the Fermi energy is highly anisotropic and its associated symmetry 
is similar to the d-wave symmetry. The anisotropy is determined by the flatness of the band dispersion. 

We emphasize that the pseudogap formation discussed is highly dependent on the specific character (perfect nesting) 
of the free band dispersion or the hopping term in the square lattice. Thus, the pseudogap is physically different from 
a gap formation such as the Mott insulator transition which takes place in a strong coupling regime. Our pseudogap 
can take place even in the weak coupling regime. 

We have focused on the undoped single-band Hubbard model with nearest neighbor hopping. Unfortunately, the 
model we have used may be quite simple to make a quantitative comparison with the experimental data. Nevertheless, 
it is possible to apply our argument in the pseudogap originated, for example, in the Hubbard model with the next 
nearest neighbor hopping (t-t'-U model). In this model, the nesting condition associated with the momentum Q holds 
only around the (tt, 0) region (and, of course, also at the other symmetric three parts in the Brillouin zone). Hence, 
the pseudogap will open in this region. However, the region in the k x = ±k y direction on the Fermi surface does 
not satisfy the nesting condition. Thus, we can predict that the electronic states around this diagonal region will 
remain in the Fermi liquid regime. Our argument is consistent with the physical picture of the so-called hot and cold 
quasiparticles by B. Stojkovic and Pines p6[ . Investigation with more realistic models will be done in future works. 
Very recently, the information of the self-energy has been extracted from the angle-resolved photoemission data 
. In the work, it has been obtained that the observed normal state pseudogap accompanies a sharp negative peak 
in ImE^ and a positive slope in ReE^. The tendency of their data is qualitatively in good agreement with our results. 

In the present work we have used a paramagnon-theory self-energy to calculate the electronic states. This treatment 
has a restriction due to the Stoner criterion in the RPA susceptibility. However, we believe that already in the level of 
this scheme, it contains important ingredients for the pseudogap formation and the partial destruction of the Fermi- 
liquid quasiparticles at the Fermi level. In particular, the anisotropic pseudogap formation in the spectral function 
in the calculated band dispersion indicates the coexistence of the Fermi-liquid-like quasiparticles and the SDW-like 
quasiparticles at the Fermi level. 
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APPENDIX: A RELATION IN THE KRAMERS-KRONIG TRANSFORMATION 

In this appendix we give a simple proof of a relation obtained from the Kramers-Kronig transformation. We apply 
this relation to the real and imaginary parts of the self-energy in the main text of the present paper. The relation 
holds for any smooth functions. Suppose the following Kramers-Kronig relation for such two functions g(x) and h(x), 

, . P f°° , , h(x') 
g(x) = - / dx'-^. (Al) 

Then, the relation can be mentioned as follows. If h(x) has a maximum (minimum) at a certain point x = xq, g(x) 
has a positive (negative) slope at xq. We assume that h.(x) has a peak at x = xq and it can be expanded around Xq 
as h(ir) ~ h(xo) + (7h/2)(a; — xq) 2 . This is reasonable as far as h(z) is a smooth function. Here, 7h is the second 
derivative of h(x) at x = xq. If 7h > (< 0), then the peak is a minimum (maximum). 

Let's observe the difference between g(xo + e) and g(xo) with e being a small shift from xq. We obtain that 

, s p f°° , / eh(a/) 

g{x + e) - g{x ) = - dx'— ; -. A2 

™ J-oo {x' - x Q y - s(x' - x ) 

We easily see that the function e/[(x' — Xq) 2 — e(x' — xq)] decreases rapidly as x' deviates from Xq and Xq + e. We 
can approximate the integration using the expansion of h(x) around x as 

g(x + e) - g(x ) ~ - dx {x , _ Xq)2 _ £{x , _ ^ ■ (A3) 

where x\ < xq < x<i, and besides, \xi — x§\, \x\ — xq\ ^> \e\. We have chosen x\ and X2 in a way which the expansion 
of h(x) is valid. By evaluating the integral up to the leading order of e, we obtain 

g(x +e) - g(x ) ^ (x 2 - xi)7h (M] 

Thus, for 7h > (7h < 0) at x — xq, the slope of g(xo) is positive (negative). 
By applying the same argument to the opposite transformation defined as 

h(a;) = -- r dx'4^, (A5) 

77 J-oo X' -X 

we find that if g(x) has a maximum (minimum) at x — xq, the slope of h(x) at x = xq is negative (positive). 

In the above derivation, we assumed that the functions of our interests can be approximated with a quadratic form 
around a peak point we have interests. One can show the same relation in the case that the function behaves as 
h(x) oc | a; — Xqj around x — xq. This is what happens if h(x) is the imaginary part of the self-energy of a marginal 
Fermi liquid [|35| . Even if the peak is a 8- function, the relation we argued is still applicable although the slope becomes 
infinite. One example of this can be seen in the mean-field Green's function for the spin-density-wave states. 
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FIG. 1. U dependence of the density of states N(u) at T/t = 0.5. 



FIG. 2. Temperature dependence of the density of states N(ui) for U/t = 2.0. Inset shows the overall structures of N{oj). 



FIG. 3. Temperature evolution of (a) the spectral function and (b) the real and (c) imaginary parts of the self-energy for 
U/t = 0.2 and at k F = (ir, 0) 



FIG. 4. Temperature dependence of (a) Inrx s (Q, v) and (b) the integrand 7k F (Q, w) of the imaginary part of the self-energy 
for U/t = 2.0. 



FIG. 5. U dependence of the excitation spectra of (a) the model susceptibility (see the text) ImxKQ,^) and (b), (c) the 
integrand i£ (Q,a;) of the imaginary part of the self-energy for the model susceptibility. 
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FIG. 6. (a) Band dispersion for several momentum directions and (b) A{k, u>) around {it, 0) region calculated for U/t = 2.0 
and at T/t = 0.22. 



FIG. 7. (a) Spectral function at Icf = {it, 0) (solid line) and (7r/2,7r/2) (dotted line), and the corresponding (b) real and 
(c) imaginary parts of the self-energy. The parameters are U/t = 2.0 and T/t = 0.22. 

FIG. 8. Spectral function at k F = {it, 0) (solid line) and {it/2, it/2) (dotted line) for U/t = 3.0 and at T/t = 0.5. Note the 
differences in scale compared with Fig. 7 (a). 

FIG. 9. Main contribution area (shown with shade) to the q summation in ImS(kF,0) for T/t = 0.22 and at (a) 
kF = (7r/2,7r/2) and (b) Isf = {it, 0) . Momenta q's in the region satisfy |/(£k F + q ) +n(£k F + q )| > 0.5. 
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Fig.l: T. Saikawa and A. Ferraz 




Fig. 2: T. Saikawa and A. Ferraz 




Fig. 3: T. Saikawa and A. Ferraz 
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Fig.4: T. Saikawa and A. Ferraz 




Fig.5: T. Saikawa and A. Ferraz 




Fig. 6: T. Saikawa and A. Ferraz 
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Fig. 7: T. Saikawa and A. Ferraz 




Fig. 8: T. Saikawa and A. Ferraz 




Fig. 9: T. Saikawa and A. Ferraz 



